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We present analytical results for all six-photon helicity amplitudes. For the computation of this loop induced process two 
recently developed methods, based on form factor decomposition and on multiple cuts, have been used. We obtain compact 
results, demonstrating the applicability of both methods to one-loop amplitudes relevant to precision collider phenomenology. 



The self-interaction of photons (light-by-light scatter- 
ing) mediated through a virtual charged fermion loop 
is a fundamental, although yet unobserved, prediction 
of quantum electrodynamics. The corresponding multi- 
photon scattering amplitudes are of outstanding theoreti- 
cal interest, since they exhibit a high degree of symmetry. 
They can be used to establish and further develop meth- 
ods for the calculation of virtual corrections to multi-leg 
processes, and to study symmetry patterns in the results, 
thus providing further insight into the analytical structure 
of quantum field theories at the loop level. 

In the past, the four-photon amplitudes were derived 
at one loop for massless and massive fermions [ [T] , and 
at two loops for massless fermions [ [5]. At two loops, 
four-photon scattering in supcrsymmctric Yang-Mills the- 
ories [[3] was studied as well, yielding first evidence for a 
leading transccndentality behaviour, which was uncovered 
subsequently also in multi-gluon amplitudes [ [4] , and has 
sparked many new developments (see [ [5] and references 
therein) in the study of multi-loop scattering amplitudes 
in super- Yang-Mills theories and perturbative gravity. 

The computation of one-loop multi-particle amplitudes 
is currently among the most pressing issues in the prepa- 
ration of precision next-to-leading order (NLO) calcula- 
tions for the upcoming CERN LHC experiments. Given 
the large variety of potentially interesting multi-particle 
final states, automated methods for one- loop corrections 
would be very desirable, and are currently under intense 
development [ EJ El HU El El QH El El El El 
EH2Dll2IlI21l23l21I23l2ai23- These methods range 
from purely analytical schemes to completely numerical 
approaches. 

Compact analytical expressions for amplitudes involv- 
ing n > 4 external photons were obtained only for specific 



helicity configurations. All amplitudes with odd n vanish 
due to parity conservation; amplitudes with even n > 4 
vanish if all or all but one photons have the same helic- 
ity [ I28| ■ For n = 6 external photons, one finds there- 
fore only two non- vanishing amplitudes, which we denote 
by A 6 ( h + ++) and A 6 ( — I 1 h). An analyt- 
ical expression for A 6 ( h + ++) was computed al- 
ready long ago [ESj, using the method described in [EE]- 

Aq( 1 1 h) was obtained recently using a purely 

numerical method for the loop integration. In this paper, 
we use two completely different recently developed meth- 
ods (based either on form factor decomposition [[6] or on 

multiple cuts [EIE1E1E]) to compute A 6 ( h + ++) 

and Aq( 1 1 (-), obtaining compact results which 

respect the symmetry properties of the process under con- 
sideration. Besides allowing a prediction for the pro- 
cess 77 — > 47, our results serve as a highly non-trivial 
proof of applicability of both methods used, and illus- 
trate how form factor-based and cut-based techniques can 
be matched onto each other in detail. 



1. Structure of the Amplitudes 

Despite the absence of a corresponding tree-level 
process and the (Feynman) diagram-by-diagram UV- 
finitcness in four dimensions, the cut-constructibility of 
the 6-photon amplitudes is not guaranteed, in accordance 
with the power counting argument in [|3J. The fact that 
the rational parts of the six-photon amplitudes actually 
do evaluate to zero was shown in [ 130] . Consequently, 
the amplitude can be written as a linear combination 
of poly-logarithms and transcendental constants, associ- 
ated to a known basis of functions, the master integrals 
(MI), formed by box-, triangle- and bubble- type integrals 
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where 

Fi = diXFi (i = l,2B). 

-Fi, -F2B are the finite parts of the one-mass and two-mass- 
easy box functions [ 1321 133] , and c?i and d 2 s are their co- 
efficients, that can be obtained from ([6]) and ([7]) below. 
The sum in Eq.(fT]) runs over the discrete quotient group 
5 4 (3,4,5,6)/(Z 2 (3,4) x Z 2 (5,6)) which has 6 elements. 
The quotient space structure can be infered from the sym- 
metry properties of the combination of basis functions. It 
is invariant under interchanging the indices er 3 «-> 04 and 
er 5 «-> erg. The permutations generate exactly those func- 
tions which are allowed by the cutting rules. 

The result for the amplitude Aq( — I 1 h) has the 

following structure 
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where 

T l = d i x Fi (i = 1, 2 A) , and Z 3 = c 3 x 7 3 . 

-Fi) Fia are the finite parts of the one-mass and two-mass- 
hard box functions, 7 3 is the triangle function with 3 off- 
shell legs and di, d 2j 4, and C3 are obtained from (Til))) . (Ti"2")) 
and given in section El The different permutation 
groups in ©, i.e. <S 3 (1, 3, 5)/Z 2 (l, 3) x 5 3 (2, 4, 6)/Z 2 (4, 6) 



for the first term, S 3 (l, 3, 5)/Z 2 (3, 5) x S 3 (2, 4, 6)/Z 2 (2, 6) 
for the second term, iS 3 (l,3, 5) x 53(2,4,6) for the third 
one and <S 3 (2, 4, 6) for the last term, contain 9, 9, 36 and 
6 elements respectively. They generate all possible differ- 
ent basis functions Fx, F2A and I3 which are allowed by 
cutting rules. Note that each cut is corresponding to a 
certain Madelstam variable. Compatibility with cutting 
rules means that no Mandelstam variables Sij or Sijk ap- 
pear as a function argument such that the corrsponding 
legs i,j, k have like-sign helicitics. This is certainly only 
true for amplitudes with massless particles. 

2. Cut-Construction 

As the six-photon amplitudes are cut-constructible, the 
computational effort is reduced to the computation of the 
rational coefficients of a linear combination of master inte- 
grals as described in section Q] According to the principle 
of unitarity-based methods [ El] , the exploitation of the 
unitarity-cuts of each master integral enables the extrac- 
tion of the corresponding coefficient from the amplitude. 
To that aim, we employ the quadruple-cut technique [ 
[14] for box-coefficients, the triple-cut integration [ [17] 
for triangle-coefficients, and the double-cut integration [ 
fT5] IT6] for bubble-coefficients, by sewing in the multiple 
cuts the QED tree-level amplitudes given in [ E3] . Spinor 
algebra and numerical evaluation of spinor products has 
been implemented in a Mathematica package [136). In the 
following we use by now standard spinor notation, with 
= ( a ~\Pi...j\b~), and multi-particle momenta 
defined as Pi...j = Pi + . ■ -+Pj, with pk being the momen- 
tum of the fc-th external photon, considered incoming. 

2.1. Construction of A 6 ( h + ++) 

The analytic expression for A§{ h + ++) was com- 
puted by Mahlon [ EH] • It can be expressed in terms of 
the two classes of box functions F\ and F^b as shown 
in eq. ([1]). Using the symmetry of the amplitude, it is 
sufficient to compute the coefficients of a representative 
box-function for each of the two classes (and their parity- 
conjugates). 
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Figure 1. Quadruple-Cuts, (left) and d^s ( rl g n t)- 
Reverse internal-helicity counterparts understood. 



displayed in Fig Q] reads 
4 X) -r!(2 3) 2 (14) 2 . 

By reversing the internal helicities one gets 
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The coefficient of Pi(s23, S34, S561), the finite part of the 
one-mass box, reads 



A , \ o (- 2 ) d T } + d i ] 
di{s2 S ,s M ,s 561 ) = 2 X — -t — L — 



(6) 



S23S34 4 

where the prefactor 2 accounts for the contribution com- 
ing from a fermion looping in the opposite direction; the 
second factor is the standard coefficient of the finite part 
of the one-mass box; and the last one is the average of 
the solutions of the quadruple-cut. The same structure is 
understood for the forthcoming four-point coefficients. 

Two-mass easy Box (6 + l"|2+|3"5+|4+) 

The coefficient of P2b(s612, S235, S6i> S35), the finite part 
of the considered two-mass-easy box is defined as 



d2s(s612, S2357 s 6l7 s 35) 

(-2) 



2 x 



d {1) + d {2) 

a 2B ' a 2B 



S612S235 ~ S61 s 35 
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with d^B being the quadruple-cut for the configuration 
displayed in FigflJ and d^ the complementary contribu- 
tion coming from a fermion circulating around the loop 
in the opposite direction. 

Although the quadruple-cuts d^g and d^s are respec- 
tively different from the quadruple-cuts an d in 
Eqs.fUE]), the coefficients of P 2 s(s6i 2 , s 235 , s 6 i, S35) in 
and of -Fi(s23, S34, S561) in © ar e not independent, and 
one finds the relation, 



d2B(s612, S235, S61, s 35 ) = -di(s 2 3, S34, s 56 i) . 



(8) 



2.2. Construction of A 6 ( — I 1 h) 

The analytic expression for A$( 1 1 h) is the 

main result of this letter. There appear three classes of 
functions: one-mass and two-mass-hard box-functions F\ 
and F2A, and three-mass triangle- function I3. Bubble- 
functions are absent. As before, we compute the co- 
efficients of a representative function for each of the 
three classes (and their parity-conjugates), and obtain the 
whole amplitude by summing over all non-identical per- 
mutations of external particles. 
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Figure 2 

di , and 
terparts understood 



From left to right: quadruple-cuts, u, 2A 




and 



di , and triple-cut c 3 ^. 



Reverse internal-helicity coun- 



Two-mass hard Box (1"2+|3"4+|5"|6 4 
We define the prefactor 



r-iA 



s 345si 6 



(26)[35](6|Pi 2 |3](2|P 3 4|5](6|P34Pi2|6)[5|P34Pi2|5] : 



(9) 



thus the result of the quadruple-cut d^l for the configu- 
ration displayed in Figgis 



d£l = r 2A s 2 U s{&l) 2 [Ab¥. 



By reversing the internal helicities one gets 
4 2 i=r 2A (l|P 34 |5] 2 <6|P 12 |4] 2 



(10) 
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The coefficient of F2a{s5&, «345i S12, S34), the finite part of 
the two-mass-hard box, reads 



^2a(s56, S345, S12, S34) = 2 



(-2) di 1 } + 4 2 i , , 
S345S56 2 



One-mass Box (1-2+3" |4+|5"|6+) 
With the prefactor 



n = 



S45 S 56S456 
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(13) 
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the result of the quadruple-cut dj 1 for the configuration 
displayed in Figgis 



4 X) =n(4 5) 2 (6|P 123 |2] 2 . 

By reversing the internal hclicities one gets 



?(2) 



r 1 (5 6) 2 (4|P 123 |2] 5 
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The coefficient of Pi(s45, S56, S123), the finite part of the 
one-mass box, reads 



di(s 45 , s 56 , s 



123) 



2 x 



(-2) 
S45S56 



(2) 
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Three-mass Triangle (1"2+|3"4 + |5"6+) 

From the triple-cut of the configuration displayed in 
Fig[2] the triangle coefficient reads 



g(i)_ 

c 3 — 

with 



(12)(l|P 34 |2] 2 e 2 N ei 



(2 4) (3 4) 



(17) 



Sl2(13)+^(si2-S34)(13) + (l|P3 4 Pl2|3)Jei 
S12S34(1 5)ei + (l|Pi2P34|5)(-Sl2 + (si2-S34)ei) 

-«ia(l 3) (2 4) + ((S12 - s 3 4)<l 3) (2 4) + 
+{34)(l|P 34 Pi2|2)-(32)(l|P34Pi 2 |4) N )e 1 l , (18) 



D„ 



s 12 (l 2) + (Js 12 -334) (1 2) + (IIP34P2 12) J ei 
si2S34(l 6)ei + (l|Pi 2 P 3 4|6)( - si2 + (si2-s 34 )ei 

-Si2<14) + ((s 12 -S34)(14) + (l|P 3 4Pl2|4))ei 
-Sl 2 (16)+((s 12 -S34)(16) + (l|P34Pl2|6))ei 
Sl2(l|P 3 4|3]ei + {l|P 12 |3]f-Sl2 + (si2-S34)ei 



- Si 2 +(si 2 -S34+(l|P34|l])ei 



where 



ei, 2 



Sl2 3S34 — S56 — S12 ± ^/ ^12,34,56 



(19) 



(20) 



2 S 12 (S 3 4-S5 6 )+S 3 4(-2S 3 4 + S 56 ) 

with 

A12, 3 4,56 = S 2 2+S3 4 + S 2 6 -2S12S34-2S12S 56 -2S3 4 S56. (21) 



The contribution coming from reversing the inner helici- 

-3 



ties, c 3 2) , 



amounts to the same value, cj? = 63 , therefore 



the coefficient of /3(si 2 , S3 4 , 85$), the three-mass triangle 
within the amplitude is, 



C3(S12, S34, S 56 ) 



2 x (4 1} 



(22) 



where the factor 2 accounts for the contribution coming 
from a fermion looping in the opposite-direction. Al- 
though not manifest in Eqs. (|18ll9p . one can see analyti- 
cally that in (fl7|) the dependence on VA drops out, be- 
cause ei and e 2 only differ in the sign of \/A. Therefore 
C3 is a rational function of spinor products. 

3. Form factor approach 

In the Fcynman-diagrammatic approach the six-photon 
amplitude is represented by 120 one-loop diagrams which 
differ only by permutations of the external photons. 

As the corresponding integrals are IR/UV finite, the 
Dirac algebra can be performed in D = 4 dimensions. 
However, using algebraic tensor reduction, one has to 
work with a (4 — 2e)-dimensional loop momentum, since 
at intermediate steps scalar integrals are generated which 
are formally divergent in D = 4 dimensions. The respec- 
tive coefficients will drop out in the end, which serves as 
a check of the computation. We use the spinor helicity 
method and define projectors on the helicity amplitudes 

Aq( h + ++) and A 6 ( 1 1 h) in such a way 

that - by choosing convenient reference momenta for the 
polarisation vectors - we obtain global spinorial factors 
for each amplitude. The resulting expressions contain in- 
tegrals with scalar products of external vectors and loop 
momenta in the numerator. For the six-point integrals 
all scalar products between the loop momentum and ex- 
ternal momenta are reducible, i.e. they can be written as 
differences of inverse propagators. As a result, at most 
rank-one six-point functions have to be evaluated. For 
the (N < 6)-point functions at most three loop momenta 
remain in the numerator. As was shown in [[30], this 
can be related to the cut-constructibility argument of [ 
[34] . To reduce the irreducible tensor integrals to scalar 
integrals we use the algebraic approach outlined in [ [5] , 
which leads to a representation of the amplitude in terms 
of the basis functions J3, Pi, P 2 a, F2B- The formalism is 
implemented using FORM [[38]. The coefficients of the ba- 
sis functions are then stored and simplified further with 
Maple and/or Mathematica. The same setup was already 
used for simpler 4- and 5-point loop amplitudes [I37j. 
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Although only a restricted set of functions actually 
needs to be evaluated due to the symmetry properties of 
the amplitude, the symmetry relations serve as a stringent 
check on the implementation. Therefore, and in order to 
test our setup in view of future applications, all 120 dia- 
grams have been calculated. 

We stress that the resulting expressions for the coeffi- 
cients, although they are not evaluated in terms of spinor 
products, allow for a fast numerical evaluation. We note 
that in all expressions at most one power of inverse Gram 
determinants survives, which is intrinsic to the chosen 
function basis. 

We have cross-checked the result obtained by the form 
factor approach with the one achieved by using the cut- 
ting rules and find perfect agreement. We have fur- 
ther compared to the recent numerical result of Nagy 
and Soper [ [23]: mapping our helicity configuration 
onto their one, according to A e (2~ , 1+, 3 — , 4 + , 6~, 5 + ) = 
Ag(l + , 2~ , 3~, 4 + , 5 + , 6~) , and using the same kinemat- 
ics as in Fig. 5 of we find the result shown in Fig. [31 
which agrees with Nagy and Soper within their plot range. 
We note that our results are produced by simply evaluat- 
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Figure 3. The modulus of the normalised six-photon am- 
plitudes s |A 6 (+H )\/a 3 and s \A 6 (-\ H — )\/a 3 

plotted for the kinematics as defined in [I23j. 



ing the analytical expressions obtained by the form factor 
approach, therefore we do not have numerical errors. The 
peak structures in the plots stem from complicated phase 
patterns in the amplitudes. 



4. Conclusions 

The six-photon amplitudes can be fully expressed as 
linear combinations of known one-loop master integrals 
with three and four external momenta. Using the tech- 
niques of form factor decomposition [ [6] and of multiple 
cuts [Q31 [111 [TBI [17], we have derived analytic expressions 
for the coefficients of the master integrals. The expres- 
sions obtained from the multiple-cut method are typically 
more compact than the form factor results, such that we 
decided to quote only the former. All coefficients agree 
numerically. We fully confirm earlier purely numerical re- 
sults [[23] for the six-photon amplitudes. Our calculation 
demonstrates the applicability of the form factor decom- 
position and multiple-cut methods to non-trivial multi-leg 
processes at one-loop, and illustrates that both methods 
can be formulated in the same integral basis. In future ap- 
plications, one could therefore envisage to combine both 
methods. Both approaches used here can also be applied 
to amplitudes containing massive particles. Our result 
shows that the analytical evaluation of one-loop ampli- 
tudes with similar kinematics relevant for the LHC is fea- 
sible. 

Note added: Shortly after the appearance of this Let- 
ter, numerical results on the same subject were presented 
by Ossola, Papadopoulos and Pittau using a reduction 
method at the integrand level [[39]. Full consistency was 
found where applicable. Moreover, the triangle coefficient 

(fT7]) in A 6 ( — | 1 h) was rederived recently by Forde [ 

[40] using an independent method. Full agreement was 
found with our result. 

Acknowledgements 

We would like to thank Jean-Philippe Guillet, Chris- 
tian Schubert and Daniel Egli for collaboration during 
early stages; Simon Badger, Ruth Britto and Andrea Fer- 
roglia, for clarifying discussions. GH would like to thank 
the ITP at Zurich University for its hospitality. This re- 
search was supported in part by Marie-Curie-EIF under 
contract MEIF-CT-2006-024178, by the Swiss National 
Science Foundation (SNF) under contract 200020-109162, 
the UK Particle Physics and Astronomy Research Council 
(PPARC) and the Scottish Universities Physics Alliance 
(SUPA). 



G 



REFERENCES 



1 



R. Karplus and M. Neumann, Phys. Rev. 80 (1950) 
380; 83 (1951) 776; B. De Tollis, Nuovo Cim. 32 
(1964) 757, 35 (1965) 1182. 

Z. Bern, A. De Freitas, L. J. Dixon, A. Ghinculov and 
H. L. Wong, JHEP 0111 (2001) 031. 
T. Binoth, E. W. N. Glover, P. Marquard and J. J. van 
der Bij, JHEP 0205 (2002) 060. 

C. Anastasiou, Z. Bern, L. J. Dixon and 

D. A. Kosower, Phys. Rev. Lett. 91 (2003) 251602. 
Z. Bern, M. Czakon, L. J. Dixon, D. A. Kosower and 
V. A. Smirnov, Phys. Rev. D 75 (2007) 085010. 
T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon and 
C. Schubert, JHEP 0510 (2005) 015; Nucl. Phys. B 
572 (2000) 361. 

C. F. Berger, Z. Bern, L. J. Dixon, D. Forde and 

D. A. Kosower, Phys. Rev. D 74 (2006) 036009. 
R. K. Ellis, W. T. Giele and G. Zanderighi, Phys. Rev. 
D 73 (2006) 014027. 

A. Denner, S. Dittmaier, Nucl. Phys. B 734 (2006) 62. 
A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, 
Nucl. Phys. B 650 (2003) 162. 

11. Y. Kuriharaand T. Kaneko, Comput. Phys. Commun. 
174 (2006) 530. 

12. C. Anastasiou and A. Daleo, JHEP 0610 (2006) 031. 

13. M. Czakon, Comput. Phys. Commun. 175 (2006) 559. 

14. R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 
725, (2005) 275. 

15. R. Britto, E. Buchbinder, F. Cachazo and B. Feng, 
Phys. Rev. D 72, (2005) 065012. 

16. R. Britto, B. Feng and P. Mastrolia, Phys. Rev. D 73, 
(2006) 105004. 

17. P. Mastrolia, Phys. Lett. B 644 (2007) 272. 

18. C. Anastasiou, R. Britto, B. Feng, Z. Kunszt 
and P. Mastrolia, Ph ys. Lett. B 645 (2007) 213; 
|arXiv:hep-ph/0612277| 



6 



7 



9. 
10 



19. R. Britto and B. Feng, arXiv:hep-ph/0612089 

20. A. Brandhuber, S. McNamara, B. J. Spence and 
G. Travaglini, JHEP 0510 (2005) 011. 

21. G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. 
Phys. B 763 (2007) 147. 

22. Z. Xiao, G. Yang and C. J. Zhu, Nucl. Phys. B 758 
(2006) 53. 

23. Z. Nagy and D. E. Soper, Phys. Rev. D 74 (2006) 
093006. 

24. R. K. Ellis, W. T. Giele and G. Zanderighi, JHEP 
0605 (2006) 027. 



25. A. Denner, S. Dittmaier, M. Roth and L. H. Wicdcrs, 
Nucl. Phys. B 724 (2005) 247. 

26. A. Lazopoulos, K. Melnikov and F. Petriello, 
|arXiv:hep-ph/0703273| 

27. C. Anastasiou, S. Beerli and A. Daleo, 



arXiv:hep-ph/0703282, 



28. G. Mahlon, Phys. Rev. D 49 (1994) 2197. 

29. G. Mahlon, in proceedings of Beyond the Standard 
Model IV, eds. J. Gunion, T. Han, J. Ohnemus, World 
Scientific (River Edge NJ, 1995). 
T. Binoth, J. P. Guillet and G. Heinrich, JHEP 0702 
(2007) 013. 

G. 't Hooft, and M. Veltman, Nucl. Phys. B 153 
(1979) 365. 

Z. Bern, L. J. Dixon and D. A. Kosower, Phys. Lett. 
B 302 (1993) 299 [Erratum-ibid. B 318 (1993) 649]; 
Nucl. Phys. B 412 (1994) 751. 

T. Binoth, J. P. Guillet, G. Heinrich and C. Schubert, 
Nucl. Phys. B 615 (2001) 385. 

Z. Bern, L. J. Dixon, D. C. Dunbar and 
D. A. Kosower, Nucl. Phys. B 425 (1994) 217; Nucl. 
Phys. B 435 (1995) 59. 

K. J. Ozeren, W. J. Stirling, JHEP 0511 (2005) 016. 
D. Maitre and P. Mastrolia, in preparation. 
T. Binoth, M. Ciccolini, N. Kauer and M. Kramer, 
JHEP 0612 (2006) 046; JHEP 0503 (2005) 065; 

JHEP 0402 (2004) 057. 

J. A. M. Vermaseren, |arXiv:math-ph/0010025| 
G. Ossola, C.G. Papadopoulos and R. Pittau, 
larXiv:0704. 12711 [hep-ph]. 
D. Forde. [arXi7:0704. 18 35 [hep-ph]. 



30 



31 



32 



33 



34 



35 
36 
37 



38 
39 

40 



